Phase randomization of three-wave interactions in capillary waves 
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We present new experimental results on the transition from the coherent-phase to the random- 
phase three-wave interactions in capillary waves under parametric excitation. Above the excitation 
threshold, coherent wave harmonics spectrally broaden. An increase in the pumping amplitude 
increases spectral widths of wave harmonics and eventually causes a strong decrease in the degree 
of the three-wave phase coupling. The results point to the modulation instability of capillary waves, 
which leads to breaking of continuous waves into ensembles of short-lived wavelets or envelope 
solitons, as the reason for the phase randomization of three-wave interactions. 

PACS numbers: 47.35.Pq, 47.27Cn, 47.52.+j 



Understanding instabilities of waves and nonlinear 
phenomena which lead to the formation of turbulence 
is important in a variety of nonlinear wave systems in 
plasma, fluids, solids, nonlinear optics, etc. Capillary 
waves, which belong to a short-wave branch of the sur- 
face waves (< 10 mm), play an important role in the sur- 
face wave physics in the ocean. They are also of interest 
from a general point of view as an example of nonlinear 
strongly dispersive waves. Capillary waves can interact 
via a three-wave interaction process which is possible due 
to the decay-type dispersion relation, u>k — / p) 1 / 2 ^/ 2 , 
where a is the surface tension coefficient and p is the 
fluid density. Such a dispersion allows three-wave res- 
onant conditions for the frequencies and for the wave 
numbers to be satisfied simultaneously: u = u>i + lu 2 and 
k = fci + k 2 . 

In the laboratory, capillary waves can be excited para- 
metrically in a vertically shaken container [H EJ |3j H], 
using electric fields [5], ultrasonic excitation [BJ, wave 
paddles [TJ, plungers |H| and other techniques. In some of 
these experiments, observations of fully developed broad- 
band spectra are reported [BJ [7J, while in others (e.g. 
[4]) discrete spectra are observed. In several experiments 
both discrete and continuous spectra are found depend- 
ing on the strength of the drive [TJ |5J IE] ■ The observed 
complex wave fields are often referred to as turbulence 
and the shapes of the spectra are compared with the 
predictions of weak turbulence theory of capillary waves 
[9j QUI • It should be noted that weak turbulence theory 
assumes that interacting waves are weakly nonlinear and 
that their phases are almost random [TJ] . However, the 
latter assumption has never been tested in experiments. 
It is also not clear if the three- wave processes lead to the 
spreading of spectral energy into broadband turbulence 
in such systems. A large body of work has also focused 
on the pattern formation in capillary waves (e.g. |12l 113] 
and Q3] for review), a topic which is not addressed here. 
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In this Letter we focus on the physics of the formation 
of turbulence driven by parametrically excited monochro- 
matic capillary waves. We show that the generation of 
continuous broadband frequency spectra occurs via non- 
linear spectral broadening of discrete harmonics. Such 
broadening results from the modulation instability of 
capillary waves, a four-wave process of sideband gener- 
ation. This nonlinear broadening changes the nature of 
the three-wave interactions between harmonics from co- 
herent to the random-phase. The phenomenology of the 
nonlinear harmonics broadening is similar to those ob- 
served during parametric excitation of the second sound 
waves in superfluid helium He 4 |15l 116] and in spin waves 
excited in ferrites [17] . 

The experiments are performed in vertically shaken 
cylindrical containers (100-200 mm in diameter, 30-50 
mm deep). Above a certain acceleration threshold capil- 
lary waves on the surface of the liquid (distilled water) are 
generated via parametric excitation of the Faraday waves 
( [HI HB] and also [Jj5] for a review) . The surface perturba- 
tions are detected using the reflection of a laser beam off 
the water surface, similarly to the method described in 
[8 . The laser beam of 5 mm diameter is reflected onto a 
diffusive screen which is imaged into a photo-multiplier 
tube (PMT) and digitized at 100 kHz sampling rate. The 
reflected light intensity is related to the tilt angle Vr? 
which is the gradient of the surface elevation rj(r,t). The 
time-varying signal of the PMT is proportional to the 
fluctuating part of the reflected beam power p(t). In the 
reported experiments, a monochromatic excitation in the 
range of the shaker frequencies of /o = 40 to 3500 Hz is 
used. This leads to the generation of the main paramet- 
rically driven f /2 subharmonic wave and a large number 
of its harmonics /„ = n/o/2. 

Nonlinear broadening of capillary waves is studied un- 
der high-frequency wave excitation at /o — 3.2 kHz. High 
excitation frequency ensures that the frequency gaps be- 
tween adjacent harmonics are substantially larger than 
the harmonic spectral widths and no overlapping occurs. 
Figure [I] shows the power spectra S(f) of the signal p(t). 
Spectra are obtained at different vertical accelerations A 




FIG. 1: (a) Fourier power spectra of the capillary waves at 
different accelerations AA — 3g, 17g, 33g, 41g, 65g. The first 
harmonic of the parametrically excited wave at /i = fo/2 = 
1.6 kHz and the zero-frequency sideband are shown. Dotted 
line corresponds to a sech 2 -fit to the exponential spectrum, 
(b) The same spectrum in the range / = (0 — 20)kHz at 
AA = 65g. (c) Spectral width Af oc 1/t as a function of the 
harmonic frequency for different accelerations. 

in the range AA = 3 to 70g, where AA = A — A t hr 
denotes acceleration above the threshold of parametric 
instability. Close to the threshold, at AA = 3g, a nar- 
row peak at f\ = 1.6 kHz is observed. This wave gen- 
erates over 10 harmonics, as seen in Figure [ljb). As 
the acceleration increases, all harmonics broaden, de- 
veloping exponential tails in the power spectrum which 
can be approximated by the squared hyperbolic secant 
function sech 2 [6(/ - /„)]. Here, S(f) = F } F* is the 
Fourier transform of p(t) and * denotes a complex con- 
jugate. As seen in Fig. [lja) exponential tails in the 
spectra extend over four orders of magnitude in £(/). 
The inverse Fourier transform of Ff = sech[6(/ — f n )] is 
s(t) = (7r/6)sech(7r 2 /& t)e l ^ nt which, in the time domain, 
describes a harmonic wave of the frequency /„ modu- 
lated by a sec/i-envelope of width r = b/n 2 . The spectral 
broadening Af cx 1/t of the wave harmonic f n increases 
approximately linear with the increase in AA. The in- 
creased input energy leads to spectral broadening, but 
does not increase the harmonic amplitudes, similar to 
other parametric systems [T^l E]- The spectral width 
also increases as a function of the harmonic number at 
fixed A A, as shown in Fig. [ljc). As the wave spectra 
broaden, a zero frequency sideband develops, Fig. [lja), 
whose width also increases with the increase in the drive. 

The development of the / = sideband indicates that 
the waves are amplitude modulated. Fig. [2] shows the 
waveforms of the band-pass filtered signals of the spec- 
trally broadened first subharmonic, / = (1.6 ± 0.2) kHz 
corresponding to the spectra of Fig. [T] The amplitude 
modulation of the waves is present even very close to the 
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FIG. 2: Waveforms of the bandpass-filtered signals (/ = 
1.6±0.2 kHz) of the first (sub)harmonic measured at different 
acceleration levels: (a,b) A A — 3g ; (c,d) A A — 33g;(e,f) 
AA — 65g. Dashed lines in (b,d,f) correspond to the sech(t) 
fit of the wavelet envelopes. 

threshold of parametric excitation. As the acceleration is 
increased, the modulation increases showing bursts of dif- 
ferent amplitudes and durations, Fig.[2ja,c,e). The shape 
of these modulation envelopes is well approximated by a 
hyperbolic secant, s(t) = sech(7r 2 /6 t). This fit is illus- 
trated as a dashed line in the right column of Fig.[2jb,d,f) 
which shows zoomed-in waveforms of the corresponding 
signals from the left column. The modulated wave ap- 
pears as a sequence of the sec/i-modulated wavelets of 
various amplitudes and widths. Thus, the frequency- 
domain and the time-domain representations are consis- 
tent with the amplitude modulation of the wave harmon- 
ics. 

The observed amplitude modulation of capillary waves 
is likely to result from the modulation instability whose 
linear theory was presented over 40 years ago [5D] but has 
not been studied in any detail, in contrast to the mod- 
ulation instability of the surface gravity waves, or the 
Benjamin- Feir instability |21j . The Lighthill criterion of 
the modulation instability, (duj/d\a\ 2 )(d 2 uik/dk 2 ) < 0, is 
satisfied for capillary waves as well as for gravity waves. 
Here, a is the wave amplitude and u> — Wfe[l — (fca) 2 /16] 
includes a nonlinear frequency correction |20j . A well- 
known exact solution of the nonlinear Schrodinger equa- 
tion which describes the evolution of the modulationally 
unstable waves jJD] is the hyperbolic secant envelope soli- 
ton, found in many physical systems including the gravity 
surface waves [22| . The above observation of the sech- 
modulated wavelets points to a similar phenomenon in 
capillary waves. 

The nonlinear spectral broadening of parametrically 
excited waves described above has also been reported in 
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FIG. 3: (a) Discrete and (b) continuous spectra of cap- 
illary waves driven by the parametrically excited wave at 
fx = fo/2 = 20 Hz. The r.m.s. acceleration is (a) AA — 2.5g 
and (b) AA = 15g. 

the second sound waves in liquid helium He 4 [15 and in 
spin waves excited in ferrites |17j . In both systems, above 
the threshold of the parametric excitation, waves develop 
exponential frequency spectra. Similarly to our observa- 
tions, it was found that the spectral width and wave am- 
plitude modulation frequency of the second sound wave 
increases with the drive |15j (the amplitude of the first 
sound wave in their case) , similar to our result of Fig. [2] 
Theoretical interpretation of the spectral broadening in 
these experiments is given in terms of the four-wave scat- 
tering [17j [23] , or a two-step three- wave interaction pro- 
cess [24] . 

Next we show results on the low frequency wave exci- 
tation. In this case, the spectral broadening of the wave 
harmonics can lead to their overlap. Figure |3ja) shows 
a V77 power spectrum of the waves excited at fo = 40 
Hz, /1 = 20 Hz at the acceleration of AA = 2.5g. More 
than n = 50 harmonics are observed in this case which 
extend over two decades in /. A spectrally broadened 
zero-frequency sideband is also seen. When the acceler- 
ation is increased to AA — lbg a continuous broadband 
spectrum is observed, Fig. (3jb). Higher acceleration, be- 
yond AA = 15<7, leads to further steepening of the surface 
ripple and to the formation of water droplets. 

As mentioned above, the dispersion relation of capil- 
lary waves allows three-wave interactions. To character- 
ize a degree of the phase coupling in discrete and con- 
tinuous spectra shown in Fig. [3] we compute the auto- 
bicoherence, or normalized bispectrum |25j . The auto- 
bispectrum of the reflected laser power p(t) is defined 

as B(/i,/a) - (FfF* fl F%) = (A^A^f-^-^), 
where / = fi + / 2 . Here angular brackets denote the 
ensemble averaging. If the phases of waves at fx, fi and 
/ are statistically random, the average value of the bis- 
pectrum is zero. The auto-bicoherence is the squared 
auto-bispectrum normalized by the auto power spectra 
of the interacting waves: 
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FIG. 4: (a,b) The auto-bicoherence and (c,d) the summed 
bicoherence of the reflected light intensity computed for the 
conditions of (a,c) discrete spectrum of Fig. |3|a), and (b,d) 
continuous spectrum of Fig. [3^" ' 



It changes between (no phase coupling) and 1 (coher- 
ent waves) and reflects the strength of the three-wave 
interactions. We also compute the summed bicoher- 
ence, SB(f) = b }(fiJ2)- This quantity gives 
f=h+h 

a measure of the total phase coupling to the frequency / 
from all frequencies fx and fi in the spectrum satisfying 
/ ./: • h- 

Figure jij shows the auto-bicoherence b 2 and the 
summed bicoherence SB (/) corresponding to the spec- 
tra in Fig. [3j The degree of the phase coupling between 
the discrete harmonics of Fig. [3ja) is very high with 
b 2 = (0.5 — 0.9) as shown in Fig. pi a). This is indica- 
tive of coherent phase interactions. In the broadband 
turbulence case of Fig. |3^b), the level of the bicoherence 
drops below 0.2 at all frequencies, Fig. |4jb), even at the 
strongest harmonics of the pumping wave. However, the 
bicoherence of the continuum increases, as seen in the cal- 
culated summed bicoherence by comparing figures Qc) 
and EJd). This suggests that at the increased excitation 
amplitude a large number of waves participate in three- 
wave interactions which have almost random phases. 

Similarly to the case of high frequency excitation, the 
wave harmonics at low frequency excitation show sech- 
spectra. Their spectral width also increases with the in- 
crease in acceleration and with the harmonic number, 
Fig.[5]ja). The increase in the spectral width results from 
the decrease in the average width r of the wavelet en- 
velopes. 

As the acceleration is gradually increased, higher har- 
monics start overlapping, such that the continuous spec- 
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ing with the increase in the drive, is linked to the de- 
velopment of the modulation instability. In the time 
domain, the instability leads to breaking of continuous 
waves into sec/i-modulated wavelets or envelope solitons. 
A decrease in the average width of the wavelet envelopes 
eventually leads to the reduction in the nonlinear inter- 
action time within resonant triads and to the drop in the 
degree of phase coupling between three waves. 

We are grateful to V. Lebedev, S. Lukaschuk, G. 
Falkovich and V. Steinberg for useful comments. This 
work was supported by the Australian Research Coun- 
cil's Discovery Projects funding scheme (DP0881571). 



FIG. 5: (a) Spectral width of the wave harmonics as a 
function of their frequency, corresponding to the spectrum 
of Fig. [5jb). (b) Power spectrum and (c) summed bicoher- 
ence of the reflected light intensity of capillary waves driven 
by the parametrically excited wave at /i = /o/2 = 20 Hz at 
intermediate acceleration of AA — Qg. 

trum is first generated at higher frequencies. This is 
shown in the spectrum of Fig. (5jb) obtained at the in- 
termediate acceleration of A^4 = Qg. In this spectrum 
the wave continuum is observed in the frequency range 
/ > 350 Hz. However the wave harmonics in the range 
/ > f c 120 Hz already lose their coherence, Fig. [5jc). 
The summed bicoherence is reduced to the random phase 
level in the range f > f c - The data of Fig. [5] (dotted 
line) suggest that the transition from coherent to random 
phases occurs at f\T < 1, where fx — 20 Hz is the fre- 
quency of the first parametric subharmonic and t « 0.04 
s is derived from the harmonic spectral width at / = f c . 
This can be explained as follows. The degree of coher- 
ence in three-wave interactions depends on the nonlinear 
interaction time between waves in the wave triads. Thus 
it is defined by the time width r of the shortest-lived 
wave (in our system the highest frequency) in the reso- 
nant triad. When this width becomes shorter than the 
period T = of the lowest frequency harmonic, at 

t/T < 1, the bicoherence is reduced. 

The experiments show that gradual development of the 
wave continuum occurs due to the spectral broadening of 
harmonics. Whether or not the Kolmogorov-type cascade 
via 3- wave interactions [51 HH] transfers energy within the 
spectrum cannot be decided solely based on the spectrum 
shape, as is often done. To further investigate this, one 
would need to analyze the wave kinetic equation using ex- 
perimentally measured wave fields, similar to the method 
described in [25] . 

Summarizing, we have shown that parametrically ex- 
cited capillary waves on the water surface become spec- 
trally broadened. The frequency spectra of the wave har- 
monics show exponential tails and are well approximated 
by the sech(6/)-fit. The increase in nonlinear broaden- 
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